..Rfsna 


NAVAL  POSTGRADUATE  SCHOX 

Monterey ,  California 


THESIS 


UTILIZATION  OF  DENSE  PACKED  PLANAR 
ACOUSTIC  ECHOSOUNDERS  TO  IDENTIFY 

TURBULENCE  STRUCTURE  IN  THE 
LOWEST  LEVELS  OF  THE  ATMOSPHERE 


by 


Louis  Robert  Moxcey 

■to;         •* 

December  1987 


Thesis  Advisor: 


D.  L  Walters 


Approved    for    public    release;     distribution     is     unlimited 


T23 


.(''^, 


J 


ECijRiTY  CLASSiF^CATiOfy  OP  TM-S   PAGE 


REPORT  DOCUMENTATION  PAGE 


a    REPORT  SECURITY  CLASSIFICATION 
JNCLASSIFIED 


lb    RESTRICTIVE   MARKINGS 


a    SECURITY  CLASSIFICATION  AUTHORITY 


b    DECLASSIFICATION 'DOWNGRADING   SCHEDULE 


3     DISTRIBUTION /AVAILABILITY  OF    REPORT 

Approved  for  public  release; 
distribution  is  unlimited 


PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


5    MONITORING  ORGANIZATION   REPORT   NUMB£R(S) 


a    NAME  OF  PERFORMING  ORGANIZATION 

4aval  Postgraduate  School 


6b    OFFICE  SYMBOL 
(If  applicable) 
61 


7a    NAME  OF  MONITORING  ORGANIZATION 

Naval  Postgraduate  School 


c    ADDRESS  (City,  State,  and  ZIP  Code) 
lonterey,    California     93943-5000 


7b    ADDRESS  (City,  State,  and  Z/PCoc/e; 
Monterey,   California     93943-5000 


a.  NAME  OF  FUNDING/ SPONSORING 
ORGANIZATION 


8b    OFFICE  SYMBOL 
(If  applicable) 


9    PROCUREMENT  INSTRUMENT  IDENTIFICATION   NUMBER 


c.  ADDRESS  (C/ty,  State,  and  ZIP  Code) 


10    SOURCE  OF  FUNDING   NUMBERS 


PROGRAM 
ELEMENT  NO. 


PROJECT 
NO 


TASK 
NO 


WORK   UNIT 
ACCESSION  NO 


1    TITLE  (Include  Security  Clarification) 

JTILIZATION   OF    DENSE    PACKED    PLANAR  ACOUSTIC    ECHOSOUNDERS    TO    IDENTIFY   TURBULENCE    STRUCTURE    IN 

:HE    lowest    levels    of   the    atmosphere     (UNCLASSIFIED) 


2    PERSONAL  AUTHOR(S) 

loxcey,    Louis   Robert 


3a    TYPE  OF  REPORT 

tasters   Thesis 


13b    TIME  COVERED 
FROM  TO 


M    DATE  OF  REPORT   (Year.  Montti,  Day) 
1987.    December 


IS    PAGE  COUNT 


73 


6    SUPPLEMENTARY  NOTATION 


COSATI  CODES 


FIELD 


GROUP 


SUB-GROUP 


19   SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 
Acoustic  Echosounders,   Atmospheric  Turbulence,    Fast  Fourier 
Transform  Algorithm 


9   ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Coherent  light  beams  propagating  through  the  atmosphere  undergo  considerable  phase 
)erturbations  due  to  fluctuating  ten^>erature  structures  in  the  atmosphere.   Understanding  emd 
leasuring  these  structures  on  a  real-time  altitude  dependent  basis  is  inherent  to  successful 
leployment  of  ground  based  l£isers  and  particle  becuns. 

One  method  used  to  detect  these  teii^)erature  structures  is  an  acoustic  profiler,  or  echo- 
lounder.   Of  immediate  interest  is  the  ability  of  high  frequency  (5  kHz)  pleuiar  array  sounders 
:o  rapidly  detect  low  level  turbulence  (below  200  meters)  and  quantify  the  results. 

This  thesis  involves  design  in^jrovements  of  previously  developed  echosounder  eurrays  aaid 
\ssociated  softwctre.   Particulary,  this  thesis  demonstrates  that  tighter  packing  of  elements 
.n  a  planau:  acoustic  azzay   produces  better  side  lobe  reduction  than  less  densely  packed  arrays 
'his  results  in  higher  energy  density  in  the  main  lobe  and  increased  performance. 

Also  included  in  this  thesis  is  a  computer  method  which  allows  relatively  accxirate  beam 
>attern  prediction  from  amy  given  plemar  array. 


20    DISTRIBUTION/ AVAILABILITY  OF  ABSTRACT 

£11  UNCLASSIFIEO'UNLIMITED       D  SAME  AS  RPT            Q  DTIC  USERS 

21     ABSTRACT  SECURITY   CLASSIFICATION                                                               1 

unclassified                                                                           | 

Ua    rjAME  OF  RESPONSBLE   INDIVIDUAL 

)onald  L.    Walters 

22b  TELEPHONE  (/nc/ude4rea  Code) 
(408)    646-2267 

^2c    OFFICE   SYMBOL                            1 

61We                                 1 

)DFORM  1473,  84  MAR 


33  APR  edition  may  be  used  until  exhausted 
All  other  editions  are  obsolete 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

0  U.S.  (Soxrnmtnt  Printlni  0"lct:   ltl(— <0('24. 


Approved  for  public  release,  distribution  is  unlimited 

Utilization  of  Dense  Packed  Planar  Acoustic  Echosounders  to 
Identify  Turbulence  Structure  in  the  Lowest  Levels  of  the  Atmosphere 

by 

Louis  Robert  Moxcey 
Lieutenant,  United  States  Navy 
B.S  ,  Jacksonville  University,   1981 

Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  PHYSICS 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
December  1987 


ABSTRACT 

Coherent    light  beams    propagating  through  the  atmosphere  undergo 
considerable  phase  perturbations  due  to    fluctuating  temperature  structures  in 
the  atmosphere.    Understanding  and  measuring  these  structures  on  a  real-time 
altitude  dependent  basis  is  inherent  to  successful  deployment  of  ground  based 
lasers  and  particle  beams. 

One  method  used  to  detect  these  temperature  structures  is  an  acoustic 
profiler,  or  echosounder.  Of  immediate  interest  is  the  ability  of  high  frequency  (5 
kHz)  planar  array  echosounders  to  rapidly  detect  low  level  turbulence  (below 
200  meters)  and  quantify  the  results. 

This  thesis  involves  design  improvements  of  previously  developed 
echosounder  arrays  and  associated  software.  Particulary,  this  thesis 
demonstrates  that  tighter  packing  of  elements  in  a  planar  acoustic  array 
produces  better  side  lobe  reduction  than  less  densely  packed  arrays.  This 
results  in  higher  energy  density  in  the  main  lobe  and  increased  performance. 

Also  included  in  this  thesis  is  a  computer  method  which  allows  relatively 
accurate  beam  pattern  prediction  from  any  given  planar  array. 
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I.  INTRODUCTION 

Optical  degradation  of  coherent  light  in  the  atmosphere  is  strongly  a 
function  of  small  scale  temperature  gradients  which  change  the  local  index  of 
refraction  in  the  atmosphere  [Ref.  1].  These  localized  structures  in  the 
atmosphere  fluctuate  rapidly,  changing  in  size,  density,  location,  velocity  and 
interaction  with  neighboring  structures.  There  is  also  considerable  variability 
with  diurnal  temperature  changes,  annual  weather  variations  and  regional 
climatological  differences.  Understanding  the  mechanisms  of  this  turbulence 
and  measuring  it  on  a  real-time  local  level  are  intrinsic  to  high  power  earth 
based  lasers. 

These  small  scale  turbulence  structures  in  the  atmosphere  play  an  even 
greater  role  with  respect  to  acoustic  absorption  and  reflection  [Ref.  2].  It  is  well 
known  that  acoustic  pulses,  on  the  order  of  1  to  6  kHz  frequency,  are  ideal  in 
defining  turbulence  structures  in  the  first  few  kilometers  of  the  atmosphere,  with 
lower  frequencies  having  greater  ranges.  Until  recently,  the  most  popular  type 
of  acoustic  sounder  was  a  monopole  speaker/receiver  mounted  in  a  parabolic 
dish  [Refs.  3  and  4].  Several  problems  with  this  method  include  difficulty  of 
transportation,  high  amplitude  side  lobes  at  large  angles  and  poor  altitude 
resolution  due  to  the  low  frequencies  involved.  These  deficiencies  can  be 
overcome  by  means  of  planar  arrays  [Refs.  5  and  6],  and  this  thesis  was 
designed  to  amplify  this  point. 


Working  under  Walters,  Weingartner  and  Wroblewski  [Rets.  7  and  8]  again 
demonstrated  the  feasibility  of  this  technique  by  building  and  testing  a  25 
element  5  kHz  planar  array.  The  result  was  an  extremely  effective  low  altitude 
echosounder  which  was  less  expensive,  more  portable  and  had  better  side 
lobe  reduction  than  previous  monople  sounders. 

During  the  course  of  this  previous  work  several  questions  arose 
concerning  the  planar  array.  If  the  elements  were  more  densely  packed,  then 
could  the  size  of  the  side  lobes  be  significantly  reduced,  thus  increasing  the 
accuracy  of  the  beam?  Could  the  array  be  modeled  on  a  computer  to  produce 
the  best  design?  Could  the  enclosure  be  improved? 

This  thesis  is  an  attempt  to  answer  these  questions,  and  involves  the 
development  of  a  dense  packed  5  kHz  hexagonal  array  and  associated  sound 
deadened  enclosure.  Additionally,  various  computer  programs  using  fast 
Fourier  transform  algorithms  were  studied  as  a  first  approximation  in 
diagnosing  acoustic  side  lobe  suppression. 

The  results  indicate  that  densely  packed  planar  acoustic  echosounders 
produce  better  side  lobe  suppression  than  more  open  arrays  and  are  accurate 
tools  in  the  identification  of  low  level  turbulence. 


II.  BACKGROUND 

The  atmosphere  induces  random  phase  perturbations  on  an 
optical  beam.  Diffraction  of  these  small  scale  phase  variations  produces  a 
larger  beam,  restricting  the  ability  to  focus  energy  on  a  target.  Optical  phase 
perturbations  in  the  atmosphere  stem  from  temperature  fluctuations  carried  by 
the  turbulent  velocity  field.  These  temperature  fluctuations  vary  in  size, 
cascading  from  tens  of  meters  in  diameter  down  to  centimeters  and  smaller. 
These  structures  appear  to  be  most  variable  and  inhomogeneous  at  the 
interface  between  stratified  layers,  particularly  near  the  earth's  surface. 

Within  the  first  few  hundred  meters  of  the  surface,  these  temperature 
structures  are  dominated  by  the  surface  heat  flux,  resulting  from  differing  rates 
of  heat  transfer  between  the  earth  and  the  atmosphere.  During  the  day  the 
surface  heats  much  more  rapidly  and  non-uniformly  than  the  atmosphere.  An 
insulating  blanket  of  air  near  the  surface  heats  and  rises  off  the  ground,  more  or 
less  in  sheets.  Warmer  sections  rise  more  rapidly,  quickly  disrupting  the 
uniformity  of  these  sheets.  Cooler  air  moves  in  to  replace  the  rising  air,  further 
disrupting  the  sheets  and  causing  vertical  thermal  plumes  to  form.  These 
thermal  plumes  continue  upwards  in  the  atmosphere,  slowly  mixing  and 
increasing  the  temperature  of  the  lower  atmosphere.  This  cycle  is  continuous, 
with  the  actual  transfer  of  energy  occurring  by  small  scale  molecular  diffusion. 
[Ref.  9]. 
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At  night  the  cycle  reverses,  with  the  surface  cooling  off  much  more  rapidly 
than  the  air  via  Planck  radiation.  In  this  case,  as  the  air  gives  up  energy  to  the 
cooler  surface,  it  produces  a  turbulent  boundary  layer  (via  velocity  interaction 
with  a  stationary  surface)  which  remains  relatively  close  to  the  ground.  This 
cooling  process  is  generally  less  turbulent  than  the  heating  process,  due  to  the 
reduced  heat  flux  at  the  boundary,  although  intense  turbulence  continues  to 
exist  at  the  interface  of  the  horizontally  stratified  layers. 

There  is  also  a  neutral  event,  during  the  transition  between  night  and  day, 
where  the  surface  temperature  equals  the  ambient  air  temperature.  During  this 
time  the  net  heat  flux  between  the  air  and  the  surface  changes  sign  and  passes 
through  a  minimum.  The  atmosphere  generally  becomes  optically  stable  during 
these  brief  pehods. 

The  structural  variability  is  clearly  evident  in  optical  instruments  designed 
to  measure  the  turbulence  of  the  atmosphere  [Refs.  10  and  11],  but  altitude 
correlation  with  optical  instruments  alone  is  difficult.  Acoustic  echosounders 
alleviate  this  difficulty. 

Acoustic  echosounders  are  devices  which  generate  and  propagate 
energy  into  the  atmosphere  and  measure  the  backscatter  return  from  turbulent 
structures  as  a  function  of  altitude.  The  power  returned  from  these  turbulent 
structures,  Pr,  is  proportional  to  the  power  transmitted,  Pj.  The  equation 
relating  these  terms  has  been  summarized  by  Neff  [Ref.  12],  based  upon  work 


by  Tatarski  [Ref  1],  and  is  repeated  in  equation  (1)  below. 

(1)  Pr=  Pj  ErEt  {ct/2)  {exp{-23R)}  (AG/R2)  h[r,/]        where, 

-  Pr  is  the  electrical  power  returned  from  range  R, 

-  Pj  is  the  electrical  power  transmitted  at  frequency  /, 

-  Er  is  the  efficiency  of  conversion  from  acoustic  to  electric  power, 

-  Ej  is  the  efficiency  of  conversion  from  electric  to  acoustic  power, 

-  c  is  the  local  speed  of  sound  (m/s), 
- 1  is  the  acoustic  pulse  length  (s), 

-  A  is  the  antenna  area  (m2) , 

-  R  is  the  range  to  the  scattering  volume  (m) 

-  exp  (-2RR)  is  the  round  trip  power  loss  due  to  attenuation  where  3  is  the 
average  attenuation  to  and  from  the  scattering  volume, 

-  G  is  the  effective  aperture  factor  for  the  antenna, 

-  |i  [r./]  is  the    scattering    cross-section  per  unit  volume;    that  is,  the 
fraction  of    incident  power  backscattered  per  unit  distance  into  unit  solid 
angle  at  frequency  /. 

The  backscattering  cross  section,  \i,  is  the  unknown  quantity  of  interest  in 
Eqn.  (1).  Tatarski  [Ref.  1]  uses  \i  in  equation  (2)  (below)  to  find  Cj^,  the 
temperature  structure  parameter. 

(2)  |i[R./]  =  0.0039  k  1/3  (Ct2  /  T2)  where, 

-  k  is  the  acoustic  wave  number  =  2k/  L, 

-  L  is  the  acoustic  wavelength, 

-  T  is  the  absolute  temperature  and 

-  Cj2  is  the  temperature  structure  parameter. 
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The  temperature  structure  parameter,  Cj^,  is  the  mean  squared 
temperature  difference  between  separate  points  in  space, 

(3)  Ct2  =  {<  T2(R2)  -  Ti  (Ri )  >2  }  /  (R2-  Ri  )2/3    ^ 

and  provides  the  scattering  centers  which  produce  backscattered  energy  from 
the  turbulence.  <AT>2  empirically  increases  by  AR2/3.  Therefore,  it  is  normalized 
by  AR2/3  to  produce  a  relatively  constant  temperature  structure  parameter. 
Since  p.  is  proportional  to  the  backscattered  return  power,  Pr,  and  to 
the  temperature  structure  parameter,  Cj2,  then  Pr  is  proportional  to  Cj2. 
Solving  equation  (1)  for  [i  and  equating  equations  (1)  and  (2),  the  following 
expression  for  Cj2  is  given  for  an  acoustic  echosounder  of  wavelength  L; 

(4)  Ct2  =  1/0.0039  (1/  ErEj)  T2  k-1/3  (2/ct)(1/AG)  Pr/Pj  {R2exp(2f3R)}. 

The  electric  to  acoustic  efficiency  factors  for  the  transducers  (speakers) 
were  measured  in  Reference  7  and  are  inherent  to  the  design  of  the  speakers. 
Er  =  Ej  =  0.5,  an  excellent  value  for  speaker  efficiency,  which  accounts  for 

their  use.   The  aperture  area.  A,  for  the  19  element  hexagonal  array  is  0.0866 
m2. 

The  effective  aperture  factor  for  the  array,  G,  accounts  for  the  non-uniform 
antenna  directivity.  G  is  a  difficult  quantity  to  calculate  or  measure.  It  represents 
a  normalized  average  of  the  acoustic  energy  within  the  transmit/receive 
antenna  beam  lobe  and  is  clearly  summarized  by  Probert-Jones  [Ref.  13]. 

1  1 


Based  upon  Reference  13  and  studies  by  Hall  and  Wescott  [Ref.  14],  a  value  of 
G  =  0.4  is  currently  used  in  computations.  This  quantity  is  presently  undergoing 
extensive  research  by  Dr.  Walters. 

Another  difficult  factor  in  Equation  (4)  involves  the  average  atmospheric 
attenuation  loss  factor,  3,  which  is  strpngiy  dependent  upon  frequency  and  the 
absolute  water  vapor  in  the  atmosphere  [Refs.  2  and  15].  This  dependence  is 
illustrated  in  Figure  1,  which  is  taken  directly  from  Reference  15.  The  water 
vapor  pressure  is  measured  via  a  digital  humidity/temperature  indicator,  model 
5165-A  from  Weathermeasure  [Appendix  A]. 
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Fig.    1.    Attenuation    versus   Water-vapor    Pressure    (mb.) 
for    Frequencies    around    1.6    KHz. 
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Combining  these  known  values  into  a  constant  Q,  equation  (4)  becomes; 
(5)  Ct2  =  (Q)  T2k-i^  {Pr  /  (ct)(  Py)}  R2  exp{23R)  where, 

Q  =  {2Pr}/{0.0039  Pt  Er  Ej  A  G}. 

Wroblewski  [Ref.  8]  utilizes  this  in  the  software  for  the  echosounderto 
calculate  a  15  minute  time  averaged  value  for  the  temperature  structure 
parameter  as  a  function  of  altitude.  Uncertainties  in  the  antenna  beam  shape, 
side  lobes,  rejection,  electronic  gains,  transducer  efficiencies,  electronic 
antenna  gains  and  the  atmospheric  attenuation  produce  uncertainties  in  the 
absolute  calibration  of  an  echosounder. 

Independent  venfication  of  the  C-pS  values  measured  and  calculated  by 
the  echosounder  are  necessary  to  confirm  the  accuracy  of  this  technique. 
Without  this  confirmation  it  is  still  possible  to  guage  relative  variations  in  the  low 
level  turbulence,  simply  with  the  echosounder  profiles.  These  profiles  are 
presented  and  discussed  in  the  following  chapters. 
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III.    EQUIPMENT 

Much  of  the  equipment  used  in  this  project  was  borrowed     from 
the    25    element    square     array     designed    by     Weingartner     and 
Wroblewski.  The  equipment  basically  involves  a  Hewlett  Packard 

(HP)  217  computer  for  data  generation  and  analysis,  an  HP  3314A 
function  generator  (controlled  by  the  software)  which  produces  the 
initial  pulse,  an  amplifier,  a  switch/pre-amp  which  passes  the 
initial  pulse  then  amplifies  and  transfers  the  return  signal  to  the 
computer,  the  array  and  shroud  which  act  as  the  physical 
transmitter/receiver  and  a  filter  used  to  eliminate  stray  noise. 
These  units  will  be  broken  down  and  individually  discussed  in  this 
section.  The  software  developed  for  the  echosounder  by  Wroblewski 
was  not  significantly  altered  and  is  not  included  in  this  document. 
[Refs.  7  and  8] 

A.       ELEMENTS 

The     elements  of  the     array     consist  of  19  piezoelectric 
Motorola  brand  speakers,     model  #  KSN  1005A,  with  attached  horns 
of     model     #  KSN   1032A.  Ceramic  piezoelectric  drivers  were  chosen 
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because  of  their  small  size,  excellent  (see  Figure  2  )  frequency 
response  at  5  kHz,  and  a  shorter  'ring'  time  with  respect  to 
electrodynamic  speakers.  Any  speaker  acts  as  a  damped 
harmonic  oscillator  with  a  'ring'  time  (recovery  time)  inversely 
proportional  to  the  damping  constant.  The  damping  is  much  stronger 
in     Piezoelectric    speakers    than    in    conventional    electrodynamic 


elements. 


Typical  Frequency  Response 
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Figure   2.      Speaker   Frequency   Response 
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An     initial    set    of    40  speakers    was    tested    to    find  19    which 
were  similarly      matched      in      frequency      response.        This     testing 
was    done    by    comparing    the     output     of     a       chosen       reference 
speaker  with      the      acoustic-to-electrical      conversion   of     each 

test  speaker.  By  clamping  the  reference  speaker  face-to-face 
with  the  test  speaker,  a  known  voltage  to  the  reference  should 
produce  a  range  of  voltages  from  the  test  specimens.  These  voltages 
were  compared  using  Lissajous  plots  on  a  Nicolet  Oscilloscope.  It 
was  more  important  to  obtain  a  set  of  identical  speakers  rather  than 
the  best  speakers,  in  order  to  achieve  a  more  uniform  final  array. 

The  horns  on  each  driver,  originally  circular  with  a  square 
flange  face  were  trimmed  to  a  circular  diameter  of  7.62cm  such  that 
they  fitted  close  packed.  Ideally  a  diameter  of  d=  L/2 
(L=wavelength)  would  provide  optimal  phase  spacing  of  elements, 
but  this  required  either  smaller  speakers  (  d  =  3.4  cm  )  which 
were  unavailable,  or  a  lower  frequency  (=2230  Hz),  which  was 
below  the  useable  range  of  available  speakers. 

Another  source  of  speaker  was  identified  as  a  Philips  brand 
piezoelectric   2002    PX     [Ref.   16]  with  better  frequency  response  at 
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3kHz  than  the  Motorola  speakers  in  use.     Unfortunately  these  Philips 
speakers    are  no  longer  manufactured. 


B.       ARRAY 


The  speaker  elements  were  configured  in  a  close-packed 
planar  array  on  an  aluminum-covered  balsa  wood  board  to  absorb 
vibrations  (Fig.  3).  The  elements  were  wired  in  parallel,  reducing 
the  array  impedance  to  300  Q  /  19  =  15.8  Q.  The  mounting  box  was 
locally  constructed  of  sheet  aluminum  and  enclosed  a  layer  of 
semi-porous  foam  loaded  with  a  thin  lead  sheet  to  reduce  stray 
noise  signals.  Input/Output  was  via  a  coaxial  connector.  Bracketts 
were  attached  to  allow  360°  rotational  testing  from  2  axes  in  the 
anechoic  chamber. 


15  cm 


FT 


Figure   3.      Array   Dimensions 
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C.       ENCLOSURE 

The  enclosure  was  designed  to  be  easily  reproducible, 
inexpensive,  portable,  and  easily  assembled  [the  earlier  enclosure 
took  some  time  to  set  up  and  break  down].  Optimally,  a  long 
cylindrical  enclosure  was  desired  to  reduce  side  lobe  scatter  and 
return  without  obstructung  the  main  lobe.  The  enclosure  chosen  to 
optimize  the  above  requirements  was  a  Rubbermaid  55  gallon  plastic 
trash  container,  with  a  height  of  84  cm,  having  a  slight  upward 
taper;  the  barrel  has  a  lower  inside  diameter  of  53  cm,  opening  to  a 
diameter  of  60  cm.  This  provided  a  finished  open  half  angle  of  about 
20°,  8°  more  than  optimally  desired  according  to  tests  in  the 
anechoic  chamber  with  the  unencumbered  array  [see  RESULTS].  A 
second  possibility  was  a  highway  barrier  barrel  [Appendix  A],  but 
these  were  much  more  expensive  and  involved  similar  dimension 
problems. 

The  inside  surface  of  the  enclosure  was  covered  with  lead  lined 
semi  porous  foam,  with  7.62  cm  convoluted  (egg  crate  shaped)  foam 
on  one  side,  and  2  cm  flat  foam  on  the  backing  side.    The  entire  foam 


18 


liner  was  attached  with  standard  contact  cement.  The  lead  lining 
amounted  to  a  density  of  about  1  Ib/sqft  which  is  sufficient  to 
reduce  one  way  noise  penetration  by  40dB.  The  convoluted  foam 
helped  break  up  internal  and  off  axis  scatter  (spectral  reflection) 
and  aided  in  side  lobe  reduction.  The  foam  used  turned  out  to  be  the 
last  of  a  salvage  batch  of  foam,  and  a  second  source  of  foam  was 
desired.  An  ideal  substitute  for  lead  lined  foam  (which  was 
difficult  to  obtain)  was  found  to  be  a  barium  impregnated  vinyl  foam 
(a  thin  rubbery  vinyl  layer  loaded  with  barium  and  sandwiched 
between  semi-porous  foam)  with  a  similar  density  of  1  Ib/sqft  and 
virtually  identical  sound  absorption  capability  (40dB).  Two 
sources  of  this  foam  were  identified  [Appendix  A]. 

D.       ELECTRONIC  HARDWARE 

The  data  processing  system  consisted  of  an  HP  Model  #217 
computer  with  a  20  megabyte  hard  drive  and  2  megabyte  memory,  an 
INFOTEK  BC  203  Basic  compiler,  analog  to  digital  (A  to  D)  converter, 
INFOTEK  FP  210  floating  point  accelerator  (to  enhance  the  speed  of 
execution),    HP    9133       floppy    disk    drive,    monitor    and    HP    ink-jet 
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printer  (see  Figure  4).  The  HP  217  was  the  workhorse  of  the  system, 
triggering  the  function  generator,  receiving  and  processing  the 
return  pulse,  and  outputting  the  results  in  graphical  mode.  Without 
this  unit,  timely  acquisition  and  reduction  of  data  would  be 
impossible. 

The  HP  3314A  function  generator,  adjusted  and  triggered 
automatically  by  the  HP  217  computer,  was  used  to  generate  a 
specific  transmitted  signal.  This  signal  was  generally  100  cycles  of 
a  sinusodial  waveform  of  1.5  volts  (rms)  amplitude. 

A  QSC  model  1700  audio  amplifier  was  used  to  increase  the 
voltage  to  the  array  from  1.5  volts  (out  of  HP  3314A)  to  30  volts 
(rms).  This  model  was  chosen  because  of  its  ability  to  drive  a 
highly  reactive  load  (the  array).  A  lighter,  more  compact  amplifier 
could  easily  replace  this  item  of  equipment,  and  will  be  considered 
during  future  upgrades. 

A  pre-amplifier  (pre-amp)  (Fig.  5)  was  used  to  isolate  the 
transmitted  pulse  from  the  received  pulse  and  to  amplify  the  return 
signal  by  a  gain  of  approximately  11,000.     The  pre-amplifer  was 
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designed  and  built  by  Walters  and  utilizes  two  low  noise  solid  state 
operational  amplifiers,  an  LT  1007  and  LT  1037.  All  four  isolation 
rectifiers  were  of  type   IN  4000. 

The  important  properties  to  consider  for  an  effective  pre-amp 
were  complete  isolation  of  the  circuit  during  transmission  (to 
prevent  overloading  components)  and  immediate  signal  retrieval 
following   the  'ring'  time  of  the  array.     The   more   rapid  the   recovery 

time  of  the  array  (a  function  of  driver  capacitance  (Cg)  and  the  loop 
resistance  (R-])),  the  quicker  the  pre-amp  can  respond  to  returned 
signals  and  the  smaller  the  surface  blind  zone.  For  this  array,  the 
capacitance  of  each  driver  was  Cq  =  0.1  jj.  F.  Thus  the  time  constant 
for  the  array  was 

RC  =  (CQ)(#drivers)(R-|)  =  time  constant 
=  (0.1*10-6  F)(19)(270  Q)  =  0.513  msec. 
For   a  given    array,    resistor   R^    provides   the   time   constant   of   the 

circuit. 

After  5  time  constants  there  is  virtually  no  'ringing'  left  on  the 
array,  and  this  sets  the  minimum  height  of  the  blind  zone.     Thus: 
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5  X  T  =  2.57  ms     and     H^jp  =   (340   m/s)(5.32*10'3s)  =  0.87  meters. 

Therefore    the  blind  zone  stands  at  approximately  1   meter  above  the 
array,  which  is  not  noticeable  on  the  echosounder  printout. 

A  Rockland  Wavetek  model  #753A  filter  {115  dB/octave}  was 
used  to  isolate  noise  and  further  amplify  the  return  signal  prior  to 
computer  analysis.  The  input  side  of  the  filter  was  set  up  as  a  low 
pass  filter  at  5  kHz  and  gain  setting  of  O  dB  (x1);  the  output  side 
was  set  up  as  a  high  pass  filter  at  5  kHz  with  a  gain  setting  of  20  dB 
(xlO).  This  produced  a  band  pass  filter  centered  at  5  kHz  with  a 
band  width  of  about  100  Hz  at  the  -3  dB  level.  The  over  all  gain  of 
the  filter  in  this  mode  was  measured  as 

output  685  m  volt 


=  9.67  =  10 


input  70.88  m  volt 

The    combination     of    the  pre-amp  and  the  filter  produced  an     overall 
return   signal   gain  of  (pre-amp  gain)(filter  gain)   =   (1 1094)(9.67) 

=  107,000. 
The  gain  of  the  system  is  thus:     107,000  *  270  a/(  270  +  15.8  Q  ) 

=    101,000. 
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IV.    RESULTS 

A.       LOBE  PATTERNS 

The  array  was  tested  in  the  anechoic  chamber  utilizing  a 
program  designed  by  Butler  [Ref.  17].  This  program  utilizes  the 
set-up  shown  in  Figure  6  and  allows  the  measurement  of  beam 
patterns  for  various  arrays.  A  band  pass  filter  was  used  to  increase 
gain  and  reduce  60  Hz  noise  from  the  equipment.  Present  on  all  the 
lobe  patterns  is  a  spurious  lobe  at  =  135°.  This  was  determined  to 
be  a  reflection  from  the  inside  door  handle  of  the  anechoic  chamber 
and  should  be  ignored.  All  lobe  patterns  have  been  normalized  to  the 
maximum  voltage  detected,  allowing  easy  comparison  between 
figures. 

The  hexagonal  array  was  rotated  around  two  axes  to  obtain  the 
maximum  sidelobe  field  of  the  array  (Figs.  7  and  8).  The  pattern  for 
the  square  array  is  shown  for  only  one  axis  of  orientation  (Fig.  9). 
Comparison  between  the  two  arrays  (Figs.  7  and  9)  shows  a  definite 
decrease  in  side  lobe  number  and  intensity  for  the  hexagonal  array 
versus  the  square  array. 
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19  ELEMENT  HEXAGONAL  ARRAY   -  APEX  MOUNTED 


MIKE 


RUN  97         MIKE-  6.2838  V       rNPUT-2. 5  V         FREQ-5050  HZ       FILTERED 


FIG.     7.         HEXAGONAL   ARRAY   BEAM  PATTERN  -  APEX  MOUNTED 
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19  ELEMENT  HEXAGONAL  ARRAY  -  SIDE  MOUNTED 


RUN  90  MIKE-  6.4951  V  INPUT-2. 5  V     FREQ-4950  HZ  355  PTS  FILTERED 


FIG.     8.         HEXAGONAL   ARRAY  BEAM  PATTERN  -  SIDE  MOUNTED 
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FIG.     9.        SQUARE    ARRAY  BEAM  PATTERN  -  SIDE  MOUNTED 
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While  the  array  was  being  tested  in  the  anechoic  chamber  a 
demonstration  of  its  detecting  ability  was  desired.  Using  a  heat 
gun,  the  array  was  set  up  in  the  chamber  exactly  as  it  would  be  in 
the  field,  except  a  Nicolet  model  #3091  digital  Oscilloscope  was 
used  to  detect  the  return  signal. 

The  heat  gun  was  tested  at  three  different  settings;  first  was 
'heat  and  air',  then  'air  only'  and  last  'heat  only'.  With  the  heat  and 
air  both  on,  a  definite  fluctuating  signal  could  be  detected  on  the 
oscilloscope  screen.  With  'air  only'  there  was  virtually  no 
detectable  signal  and  with  'heat  only*  the  signal  was  almost  back  to 
its  level  of  'heat  and  air'.  No  measurements  were  taken,  but  the 
demonstration  was  repeated  several  times.  This  directly  illustrated 
that  the  turbulence  structures  in  the  atmosphere  are  mainly  a  result 
of  temperature  fluctuations  and  not  due  to  turbulent  velocity 
fluctuations. 

The  next  step  involved  testing  the  side  lobe  suppression  of  the 
enshrouded  array.  The  mechanics  of  this  operation  were 
considerably  more  involved  than  simply  suspending  the  array  from  a 
rotating  motor.     Being  too  heavy  to  suspend,  the  shroud  was  taped  to 
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a   rotating   turntable     and   centered   beneath   the   rotator   motor.      The 
motor  was  connected  to  the  shroud  via  aluminum  rods  and  clamps. 
When    the    rotating  motor    turned,     so    did    the    shroud,  allowing  the 
enshrouded  beam  pattern  to  be  measured  and  recorded  with  the  same 
software  as  used  previously. 

The  initial  beam  pattern,  (Fig.  10)  showed  that  the  shroud,  while 
significantly   reducing  side  lobes,  was  not  long  enough  to  completely 
eliminate    them.     Some  of  the  side  lobes  were  getting  through. 

A  first-cut  step  in  further  reducing  these  side  lobes  involved 
placing  a  15  cm  wide  sleeve  of  barium  impregnated  smooth  foam 
around  the  upper  inside  circumference  of  the  shroud.  The  result  (Fig. 
11)  was  not  promising,  and  demonstrated  that  further  refinement 
was  required.  One  thing  Figure  11  demonstrated,  however,  was  that 
specular  reflection  from  a  smooth  surface  increases  side  lobe 
activity  and  thus  convoluted  foam  shrouds  (vice  smooth  foam 
shrouds)  are  required. 

The  next  step  involved  placing  a  15  cm  wide  annular  ring  coated 
with  sound  deadening  foam  over  the  exit  of  the  shroud  (similar  to  a 
trash  can  lid  with  the  center  removed).     This  left  an  open  circle  of 
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ENCLOSED  -  HEXAGONAL  ARRAY  -  FILTERED 


RUN   115       MIC  =  .157885  V       INPUT-2. 5  V         FREQ=5000  HZ 


FIG   10  ENSHROUDED  ARRAY 
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ENCLOSED  -  HEXAGONAL  ARRAY  -  FILTERED 


MIC 


RUN  MODI     MIC  -.11 4418  V       INPUT-2.  5  V         FREQ=5000  HZ 


FIG.   11.        ENSHROUDED   ARRAV  WITH  RING  LINER 
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48.25  cm  diameter  for  the  acoustic  waves  to  propagate.  This 
adjustment  slightly  decreased  the  side  lobes  (=  3  dB),  but  not 
significantly  (Fig.  12).  Further  narrowing  of  the  neck  via  more 
annular  rings  only  succeded  in  reducing  the  main  lobe,  with 
little  subsequent  change  in   side  lobe  activity. 


B.       FIELD  MEASUREMENTS 

Field  measurements  were  conducted  on  several  occassions 
over  a  six  month  period  at  various  locations  around  California.  They 
generally  coincided  with  optical  measurements  being  conducted  by 
Walters  [Refs.  9  and  10]  and  were  used  to  supplement  the  optical 
data.  Of  greatest  interest  were  mountaintops  above  the  marine 
boundary  layer  (=4000  ft.),  where  the  turbulence  should  be  a 
minimum. 

Two  mountaintops  in  particular  were  visited,  one  being  Lick 
Observatory  on  Mt.  Hamilton  and  the  other  being  Mt.  Wilson 
Observatory  at  Mt.  Wilson.  These  were  chosen  because  of  the 
availability  of  telescope  domes,  which  were  required  to  reduce  wind 
vibration  on  the  optical  components  and  protect  the  computers  from 
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FIG.  12.   ENSHROUDED  ARRAY  WITH  ANNULAR  MUZZLE 
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dust  and  debris.     The  echosounder  itself  is  impervious  to  wind,  and 
can   be  operated  wherever  electricity  is  available. 

Lick  Observatory  was  visited  on  1  and  2  August  1987  and  some 
of  the  results  are  presented  in  the  following  figures.  All  times  are 
given  as  universal  time.  Subtract  seven  hours  to  obtain  local 
(daylight  savings)  time  at  the  site.  Both  days  were  similar  with 
respect  to  turbulence.  Therefore,  only  profiles  from  2  August  will 
be  presented  to  illustrate  the  ability  of  the  echosounder  to  detect 
and  present  the  rapidly  varying  nature  of  the  atmosphere  accurately. 
The  Cy^  profiles  are  currently  unverified  and  should  be  ignored. 

Figure  13  was  taken  early  in  the  evening  and  demonstrates  the 
presence  of  convective  pluming.  This  pluming  vanishes  gradually, 
eventually  yielding  to  the  nighttime  structure  which  starts  in  Figure 
14.  By  this  time  a  definite  layer  has  set  in  and  surface  pluming  is 
no  longer  dominant.  Later  in  the  evening  a  strong  turbulence  layer 
descendes  (Fig.  15)  which  quickly  settles  [Fig.  16]  and  is  quite  dense. 
This  layer  remains  close  to  the  ground  the  rest  of  the  evening,  with 
light  turbulence  detectable  above  it. 
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Fig.  14.  Echosounder  Profile;  2000-2015  Local  Time 
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Early  in  the  morning  (Fig.  17)  another  strong  layer  is  noted.  As 
the  sun  begins  to  rise  (Fig.  18),  commencement  of  the  daily 
convective  pluming  becomes  clearly  evident.  These  figures 
illustrate  the  dramatic  ability  of  the  high  frequency  echosounder  to 
profile  low  level  turbulence  in  the  atmosphere.  Interestingly, 
whenever  a  turbulence  structure  was  detected  on  the  profiler,  there 
was  also  noticeable  degradation  to  both  the  optical  instruments  in 
use. 

The  next  occasion  to  use  the  echosounder  occurred  on  6  and  7 
November  1987  with  a  field  trip  to  Mt.  Wilson  Observatory. 
Unfortunatley,  no  echosounder  data  could  be  obtained  due  to  severe 
radio  frequency  [RF]  interference  from  approximately  10  radio/TV 
transmitters  co-located  at  the  observatory.  The  trip  was  not  a  loss, 
however,  as  it  was  used  to  correct  a  subtle  flaw  in  the  design  of  the 
pre-amp. 

The    problem  was    traced  to  the  'floating  point  ground'  design 
of  the  pre-amp.      The  coaxial  cables  leading  to  the   pre-amp   made 
ideal    antennas    at   the    200    megahertz    range,    transferring    enough 
radio  frequency  voltage  into  the  pre-amp  casing  to  influence  the 
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circuit.  The  rectified  voltage  effectively  overloaded  the  first  stage 
amplifier,  overriding  the  weak  return  signal  from  the  array  and 
destroying  the  data.  This  effect  went  undetected  in  previous  use 
of  the  pre-amp.  The     problem     was  corrected  by  improving  the 

grounding  structure  of  the  device. 

Note  that  all  these  profiles  were  taken  with  the  hexagonal 
array  mounted  in  the  square  shroud,  as  the  hexagonal  shroud  was  not 
yet  completed.  To  ensure  that  the  quality  of  data  was  not 
compromised,  the  square  array  was  also  run  at  the  site  and  the 
profiles    were    identical. 
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V.   ANTENNA  BEAM  PATTERN   MODEL 

This  thesis  also  modeled     the      planar    acoustic     array  as     an 
optical    aperture,    obtaining    the    Fourier    transform    of    the    far   field 
'diffraction'   pattern   via  computer  calculations.      The  technique   rests 
on  the  fact  that 

".  .  .  the  Fraunhofer  diffraction  pattern  is.  the  Fourier  transform 
of  the  field  distribution  across  the  aperture  .  .  .".      [Ref.  18,  p.  412] 
The  quote  applies  to  electromagnetic  radiation,   but  Dr.  waiters 
correctly   construed    the    electromagnetic   amplitude    as    an    acoustic 
pressure,  obtaining  the  following  equation; 

(6)         P{k,,ky)  =  1/iR  llgperture  ^^^'V)  ®^P  ^'  (^^x^  +  ^V))  dxdy  ; 

-  P(kx,ky)  is  the  pressure  distribution  in  the  image  plane, 

-  k^  and  ky  are  the  spatial  frequencies  along  the  x  and  y  axes  of 

the  aperture  respectively;    k|=  k*J/R,  where  J=  x  or  y;  k=27r/X  , 

-  P(x,y)  is  the  pressure  amplitude  at  the  aperture, 

-  exp  {i  (k^x    kyy)}  is  the  spatial  variation  of  the  aperture  and 

-  R  is  the  distance  to  the  image  plane;  i  accounts  for  temporal 

phase    distributions  across  the  aperture. 
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The  model  can  be  thought  of  as  a  plane  coherent  acoustic  wave 
passing  through  the  aperture,  which  produces  a  diffraction  pattern 
(the  Fourier  transform)  in  the  far  field  (the  image  plane)  by  means 
of  the  Huygens-Fresnel  principle.  From  this  pattern  it  was  possible 
to  reasonably  predict  the  size  of  the  main  lobe  and  the  approximate 
angles  and  intensities  of  the  side  lobes.  This  modeling  was  done  for 
both  the  square  array  and  the  more  densely  packed  hexagonal  array. 
The  results  are  presented  in  this  chapter. 

A.       HARDWARE 

The  computer  used  in  this  phase  was  a  Hewlett  Packard  (HP) 
217  with  an  INFOTEK  Basic  compiler,  identical  to  the  computer  used 
in  the  echosounder  arrangement.  The  memory  requirements  for  this 
program,  however,  were  tremendous,  requiring  over  4  megabytes  of 
memory  when  using  the  512  x  512  grids.  This  lead  to  insertion  of  2 
additional  memory  boards  into  the  computer,  (one  from  INFOTEK  and 
one  from  EVENTIDE),  bringing  the  hard  memory  to  over  5  megabytes. 


46 


B.       SOFTWARE 

The  program  [Appendix  B]  utilizes  a  Fast  Fourier  Transform 
(FFT)  routine  supplied  by  INFOTEK.  The  idea  is  to  enter  the  aperture 
(array)  as  a  small  series  of  points  on  a  huge  grid.  Typically  the  grid 
size  would  be  256  x  256  points  square,  and  occassionally  512  x  512. 
The  array  takes  up  only  a  small  center  portion  of  this  grid, 
approximately  40  x  40  grid  points.  The  larger  the  grid  field  in 
relation  to  the  input  aperture,  the  better  the  resolution,  provided  the 
input  aperture  is  not  so  small  as  to  obscure  the  finer  details. 

The  aperture  area  was  entered  digitally  as  a  series  of  'ones', 
signifying  a  uniform  coherent  plane  wave  at  the  aperture,  with  the 
remainder  of  the  field  entered  as  zeroes.  The  FFT  routine  then 
computed  the  far  field  Fourier  transform  of  the  array,  and  plotted 
the  output  intensity  in  relation  to  a  gray  scale.  The  gray  scale  was 
adjusted  to  encompass  3  decades  of  intensity. 


C.       RESULTS  AND  ANALYSIS 

A       sample     of     the     input       apertures,     both     square     and 
hexagonal,      are  shown  in  Figures  19  and  20.      Generally,  the  radius 
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FIG.  19.   HEXAGONAL  ARRAY  INPUT 


FIG.  :0.   SQUARE  ARRAY  INPUT 
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of  each  element  is  only  4  grid  units,  with  a  background  field  of  256 
X  256  grid  units.  These  figures  have  been  enlarged  to  show  detail. 
The  output  Fast  Fourier  Transform  for  each  is  shown  in  Figures  21 
and  22  respectively.  The  center  portion  of  the  FFT  for  each  array  is 
the  lowest  spatial  frequency  for  the  array,  and  corresponds  to 
the  main  lobe.  Rings  further  out  correspond  to  the  side  lobes  at 
larger  and  larger  angles.  Rings  far  from  the  main  lobe  are  not 
physical  except  at  very  high  frequencies.  Interestingly,  as  the 
frequency  is  increased,  more  and  more  side  lobes  are  forced  into  the 
forward    hemisphere. 

The  program  calculates  the  beam  structure  of  a  single  piston 
driver  on  an  infinite  plane.  The  actual  source  is  not  a  piston  driver 
(but  the  approximation  is  close).  Huygens  wavelettes  propagating 
backwards  from  the  end  of  the  horns  reflect  from  the  support  plane, 
increasing  the  energy  of  the  actual  lobe  field.  Also,  the  physical 
side  lobes  can  extend  no  further  than  90°  from  the  axis  of  the  array 

(due    to    surface    absorption    and    reflection)    but   this    limitation    is 
ignored    by    the    program.    In    addition,  an    obliquity    factor    was    not 
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included.  These  limitations  affect  the  graphical  output,  but  do  not 
significantly    affect    the    following    calculations. 

The  angles  of  the  output  beam  pattern  rings  depend  on  the 
frequency  of  interest  (5  kHz),  the  dimensions  of  the  actual  array  and 
the    input/output    dimension    of    the    aperture.  To       calculate 

information  from  these  printouts  it  is  necessary  to  use  optical 
aperture  equations  from  Hecht  and  Zajac  [Ref.  18].  The  half  angle 
formula     for  a  circular  aperture  is: 

(7)  sin  0  =  q/R  =  1 .22  L72a  where, 

-  q  is  the  distance  from  axis  to  side  lobe  of  interest, 

-  a  is  the  radius  of  aperture, 

-  L  is  the  wavelength, 

-  R  is  the  distance  to  the  far  field, 

-  0  is  the  half  angle  or  cone  angle  from  the  axis  to  q. 

The  hexagonal  array  is  not  quite  circular,  but  in  the  limit  it 
approaches  a  circle.  In  this  case  the  constant  1.22  is  simply 
replaced  by  F,  a  constant  to  be  determined  by  experiment. 

By  using  Brigham's  text  [Ref.  19]  it  is  possible  to  convert 
equation  (7)  into  the  form: 
(8)  0  =  F  ARCSIN   (Ln/NAT)  where. 
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-  0  is  the  angular  measure  of  the  side  lobe  from  the  main  lobe, 

-  L  is  the  wavelength  of  interest  =  6.86  cm, 

-  NAT  is  the  total  width  of  the  sampled  grid  =  (ND/X), 

-  N  is  the  width  of  the  input  window  (grid  units;1 28,256  or  512), 

-  D  is  the  actual  dimension  of  the  physical  array  (cm), 

-  X  is  the  dimension  of  the  input  array  in  grid  units, 

-  n  is  the  output  number  of  grid  units  to  the  side  lobe  of 

interest=(rz/w), 

-  r  is  the  distance  to  the  output  side  lobe  of  interest  (cm), 
{NOTE:    X  and  D  must  correspond;  r  and  n  must  correspond}, 

-  z  is  the  width  of  the  output  window  in  grid  units  (<  N/2), 

-  w  is  the  width  of  the  input  and  output  windows  in  centimeters. 

As  an   example,  consider   Figure  23.        For  this  example,    D   = 

38.735  cm   (measured     between  the     outer  edges  of  the  farthest  2 

speakers),    w  =  10.3  cm,  N  =  256  grid  points  and  z  =  128  grid  points. 

To  calculate  the  angular  extent  of  the  main  lobe  in  real  space,  first 

measure  the  distance  (in  cm)  from  the  corner  of  the  output  quadrant 

to   the   edge  of  the   main   lobe,   along   the   horizontal   axis   (Fig.   23). 

Here,   r  =  0.7  cm  and   X  =  40  grid  points.     Therefore,   equation   (8) 

becomes: 

(6.86)(0.7)(1 28/1 0.3)  cm 

(9)  0  =     F  X  ARCSIN  =  F  X  13.9°  . 

(256)(38.735)/40  cm 
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Figure  23.     Output  Quadrant  of  Hexagonal  FFT 


The  half-angle  of  the  actual  array  is  only  12°  as  measured  on 

Figure  8.    This  means  that  F  =  0.86  or  1/F  =1.16  which  is  very  close 

to    1.22.        Using   the  value   F  =   1.16,   then   the   90°   point  can   be 

calculated  for  the  output  field.     In  this  case, 

(256)(38.735) 

(10)  Sin  {(1.16)(90°)}         =     r     =     2.81   cm, 

(6.86)(128/10.3)(40) 
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which  means  everything  outside  the  corner  radius  of  2.81   cm  can  be 
ignoredas  it  is  non    physical. 

The  output  FFT's  were  correlated  to  the  actual  lobe  patterns  of 
Figures  7  and  8,  the  apex  and  side  mounted  patterns  for  the 
hexagonal  array,  and  the  results  are  plotted  in  Figures  24  and  25 
and  tabulated  in  Table  1.  The  amplitudes  were  normalized  by 
assigning  a  value  of  unity  to  the  main  lobe;  the  angles  are  measured 
counter-clockwise  from  the  standard  X  axis;  the  amplitude  scale  is 
in  decibel   units. 

The  measured  amplitude  was  calculated  using  INVLOG  (dB/20), 
where  'dB'  is  the  normalized  decibel  level  estimated  from  Figures  7 
and  8.  The  computed  amplitudes  relied  on  the  formula  20  LOG 
"Q/Max"  [see  Appendix  B,  lines  1390  and  2380  for  definitions  of  Q 
and  Max]. 

The  angles  were  computed  using  a  simpler  form  of  equation  (8); 

(11)  0    =    ArcSin  {  (K-1)  LVNAt  },  where  K  is  an  index  from 

1    to    N/2,    and   essentially   allows   the   computation    of   the   angle    at 

each   grid   point  chosen.      The   measured   angles  were  taken   directly 

from  Figures  7  and  8. 
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Fig.    24    Polar  Plot  of  Apex  Mounted  Hexagonal  Array 
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Fig.  25.  Polar  Plot  of  Side  Mounted  Hexagonal  Array 
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TABLE   1.  COMPARISONS  BETWEEN  ACTUAL  LOBE  PATTERNS  AND 
COMPUTED  LOBE  PATTERNS. 


APEX   MOUNT   (  Fig.  7  vs.  24  ) 


Measu 

ired 

Com 

puted 

Q 

amolitude 

0 

amolitude 

1  St  m  i  n 

1 1° 

.025 

12° 

.002 

1  St  max 

16° 

.126 

17° 

.128 

2nd  min 

21° 

.056 

23° 

.007 

2^^  max 

26° 

.095 

28° 

.108 

3''d  min 

33° 

>.01 

38° 

.006 

3''d  max 

42° 

.056 

41° 

.027 

4th   min 

51° 

.018 

46° 

.007 

4th  max 

64° 

.056 

62° 

.148 

SIDE   MOUNT      (  Fig.  8  vs.  25) 


Measured 

Computed 

Q 

^mplitud? 

0 

gmplitucj? 

1  St  min 

12° 

.050 

14° 

.048 

1  St  max 

18° 

.158 

19° 

.151 

2nd  min 

27° 

>.01 

28° 

.01 

2nd  max 

35° 

.039 

36° 

.030 

3''d  min 

52° 

>.01 

46° 

.024 

3''d  max 

61° 

.025 

56° 

.029 

56 


The  side  mounted  section  corresponded  to  measurements  along 
the  X  axis,  and  the  apex  to  measurements  along  the  Y  axis.  The  error 
of  the  measurements  was  grossly  estimated  as  approximately  10% 
for  the  minimums,  as  they  tended  to  be  obscured,  and  3%  or  less  for 
the  maximums.  Some  values  exceed  this  amount  considerably  and 
require    further    investigation. 

Generally,  however,  the  tabulated  values  agree  strongly  and 
suggest  that  this  technique  has  considerable  credibility  at 
predicting   side  lobe  structures  of  untested  arrays. 
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VI.    CONCLUSIONS 

This  thesis  demonstrated  that  increasing  acoustic  planar  array 
speaker  density  from  25  speakers/1452  cm^  (1.72  x  10-^ 
speakers/cm^)  to  19  speakers/943  cm^  (2.02  x  10'^  speakers/cm^) 
effectively  suppressed  the  lobe  spreading  from  a  70°  cone  angle  to 
less  than  a  45°  cone  angle,  with  a  reduction  in  side  lobe  intensity  of 
between  4  and  10  dB. 

Graphical  profiles  demonstrated  the  great  variability  and 
volatility  of  turbulence  structures  in  the  atmosphere  and  underline 
the  need  for  better  tools  to  measure  and  understand  this  structure. 
These  profiles  identify  gravity  waves,  thermal  plumes,  boundary 
layers  and  other  temperature  induced  structures.  Particulary 
important,  this  device  allowed  tentative  quantification  of  altitude, 
thickness  and  relative  intensity  of  mountain  boundary  layers  over  a 
potential  DOD  laser  site. 

Independent  verification  of  the  C-i-2  values  is  important  and  is 
currently  under  study.  It  will  further  enhance  the  utility  of  this 
device.       Due    to    the    strong    exponential    dependence    of    acoustic 
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attenuation  (3)  with  increasing  frequency,  future  improvements 
should  also  include  lower  frequency  drivers  (2-3  kHz)  of  similarly 
high  efficiency,  to  improve  the  range  of  the  device. 

It  was  shown  that  acoustic  array  beam  patterns  could  be 
successfully  modeled  on  a  computer  by  means  of  a  Fast  Fourier 
Transform  algorithm  and  tested  for  optimal  design.  Future 
improvements  to  this  technique  include  a  method  to  stochastically 
simulate  turbulence  and  model  the  beam  structure  as  it  passes 
through  this  artificial  turbulence.  Calculation  of  main  lobe  energy 
versus  total  energy  for  various  arrays  should  also  be  studied. 
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APPENDIX     A 


List    of   Material    Sources 


Disclaimer:  These       sources    are    listed    only    in    case       exact 

duplication  of  this  thesis  is  required  at  some  future  date.  Other 
equally  competitive  sources  for  each  item  do  exist  and  their 
exclusion  is  not  a  reflection  of  product  quality  in  any  way. 


Barium  impregnated  foam: 
1.    Barrier  Corp. 
Tigard,  Oregon 
(503)    639-4192 
ROC:  Mark  Dove 


2.    Latchford  Corp. 
(714)  734-3660 
ROC:   Tom  Relligreeni 


Cylindrical    Enclosure: 
1.    Rubbermaid,  Inc. 
1147  Akron  Rd 
Wooster,  Ohio    44691 
(216)    264-6464 


2.     Traffic  Products  Co. 
Citrus  Heights,  Ca 
(916)  723-2958 
ROC:   Michael  Fish 


Speakers: 

Motorola  Inc. 
4800  Alameda  Blvd  N.E. 
Albuquerque,  N.  M.    87113 
(505)    822-8801 

Temperature/humidity    Indicator 
Weathermeasure 
P.O.  Box  41039 
Sacramento,  CA    95841 
(916)    481-7565 
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APPENDIX     B 
Fast  Fourier  Transform  Program 


10  !  Progran  plots  either  the  input  p-oints,  the  FFT  or  the 

Polar  Plots  of  a  19  element  heAagon  array  or  a  Z5 

20  !  element  square  array . 

30  !   This  program  is  interactive. 

40  I  dlbj  Dec.  18,  8G/LRM  modified  in  OCT  87 

50  COM  Re(  t  :512  )  ,Im(  1  :512  ) 

S0  INTEGER  M,N 

70  DIM  Real(512  ,512  )  ,Imag(5l2 ,512  ) 

80  REnL  Amp(  100  )  ,Ang(  100  ) 

90  1 

100  PRINT  "ENTER  THE  GRID  SIZE  N;  N=  512,256  OR  128" 

1 10  INPUT  N 

120  PRINT 

150  I 

140  N2-N/2 

150  M-L0G(N)/L06(2  ) 

180  I 

170  PRINT  "ENTER  ' 1 '  TO  SEE  THE  19  ELEMENT  HEXAGON." 

180  PRINT  "ENTER  '2'  TO  SEE  THE  25  ELEMENT  SQUARE." 

190  IfJPUT  X 

200  PRINT 

210  PRINT 

220  PRINT  "ENTER  '1'  TO  SEE  A  PLOT  OF  THE  INPUT  POINTS" 

230  PRINT  "ENTER  '2'  TO  SEE  A  PLOT  OF  THE  FFT" 

2^0  PRINT  "ENTER  '3'  TO  SEE  A  POLAR  PLOT  OF  INTENSITY" 

250  INPUT  A 

2G0  PRINT 

270  PRINT 

280  PRINT  "ENTER  THE  GRID  SIZE  OF  EACH  ELEMENT,  R." 

290  PRINT  "  USUALLY  R=4  OR  8." 

300  INPUT  R 

310  PRINT 

320  PRINT 

330  PRINT  "ENTER  '7'  TO  SEE  THE  FIRST  QUADRANT  ONLY" 

340  PRINT  "  ELSE  ENTER  '0' .  " 

350  INPUT  C 

3S0  I  B=  THE  SIZE  OF  THE  OUTPUT  UINDOU 

370  B=N2/2 

380  ! 

390  IF  X=l  THEN 

400  I 

410  F  =  2.  1 ♦R*C05(  .5235988  ) 

420  FOR    I=N2-6«R    TO    N2+6»R 
430  FOR    J=N2-6*R    TO    N2+6»R 

440  IF    (  (  I+4»R-N2  )"2  +  (  J-N2  )"2  )-''R*R    THEN 

450  ReaK  I  ,J  )=1  .000  ■ 


460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
G50 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 
830 
890 
900 
910 
920 
930 
940 
950 
960 


FOR  K=1  TO  4 

ReaK  I+2»K*R,J  )=Real<  I  ,J  ) 
ReaK  I+2»K*R-R  ,J  +  F  )  =  Real(  I  ,J  ) 
ReaK  I  +  2*K*R-R  ,J-F  )  =  Reai(  I.J) 
NEXT  K 
FOR  K=l  TO  3 

ReaK  I  +  2*K*R,J-2*F  )=Real(  I  ,J  ) 
ReaK  1+2*  K*R,J  +  2»F  )=ReaK  I  ,J  ) 
NEXT  K 
END  IF 
NEXT  J 
NEXT  I 
I 

EfJD  IF 
I 

IF  X=2  THEN 
I 

K-N2+5»R 
L=N2-5*R 
FOR  1=1  TO  N 
FOR  J=l  TO  N 

IF  KK  AND  I>L  AND  J'K  AND  J>L  THEN 
Cy=N2-4*R 
FOR  P=l  TO  5 
Cx=N2-4*R 
FOR  0=1  TO  5 

IF  IvCx  +  R  AND  KC^.-R  AND  J^-Cy-R  AND  JxCy  +  R  THEN 
IF  R»R/(  (  I-Cx  )-^2  +  (  J-Cy  )'^2  )  THEN 

ReaK  J  ,1  )  =  1  .000 
END  IF 
END  IF 
Cx=Ca+2*R 
NEXT  0 
Cy=Cy+2*R 
NEXT  P 
■  END  IF 
NEXT  J 
NEXT  I 
I 

END  IF 

PRINTER  IS  701 
I 

!  This  section  plots  the  input  points  only 
I 

IF  A=l  THEN  GOTO  2240 
I 

I   THIS  SECTION  IS  THE  FAST  FOURIER  TRANSFORri 
( 

TI=TIMEDATE 
FOR  1=1  TO  N 
FOR  J=l  TO  N 
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970        Re( J  )  =  Real(  I  ,J  ) 

980        In( J  )  =  Imag(  I  ,J  ) 

990      NEXT  J 

1000  I 

1010     CALL  Ff t(M  ) 

1020  I 

1030     FOR  J=l  TO  N 

1040      ReaK  I  ,J  )  =  Re(  J  ) 

1050       Imdg( I  ,J  )  =  In( J  ) 

10S0     NEXT  J 

1070  NEXT  I 

1080  < 

10^0  FOR  J=l  TO  N 

I  100     FOR  I-I  TO  N 

1110      Re( I  )=Real(  I  .J  ) 

1  120       Im(  I  )  =  Imag(  I  ,J  ) 

1 130     NEXT  I 

1140  i 

1  150     ChLL  Ff  t(ri  ) 

1  1G0  I 

1 170     FOR  1=1  TO  N 

1  180      ReaK  I  ,J  )  =  Re(  I  ) 

1  1  90       Imdg(  I  ,J  >-Ifvi(  I  ) 

1200     NEXT  I 

1210  NEXT  J 

1220  T2-TINEDATE 

1230  I  PRINT  INT(  100(T2-T1  )/60  )/100, "MINUTES" 

1240  I 

1250  IF  A=2  THEN  GOTO  1960 

1260  I 

1270  I  This  section  graphs  the  polar  plot  of  the  FFT. 

1280  I 

1290  Bl=N2/2 

1300  B2=1 

1310  '  D=  The  actual  dimension  of  the  array(cpi) 

1320  0=38.735 

1330  PnaA  =  (  1+INT(D*N/(68.5»R  )  )  ) 

I  340  I 

1350  '  The  'K'  loop  does  the  side  and  ape.-,  mounts 

1360  I 

1370  Rl=Real( 1,1) 

1 380  Ia=Imag( 1,1) 

1390  MaA  =  SQR(Rl*Rl  +  Ia*Ia  ) 

1400  I 

1410  FOR  K=l  TO  2 

1420  I 

1430  IF  X=l  AND  K=l  THEN 

1440     PRINT  '■  HEXAGONAL  ARRAY  ,  SIDE  flOUNTED" 

1450  END  IF 

1460  IF  X=l  AND  K=2  THEN 

1470     PRINT  "  HEXAGONAL  ARRAY,  APEX  MOUfJTED" 
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14  80  END  IF 

1490  IF  X=2  AND  K=l  THEN 

1500     PRINT  "  SQUARE  ARRAY,  SIDE  MOUNTED" 

1510  END  IF 

1520  ! 

1530  '  This  section  computes  the  ape^  square  array 

1540  I 

1550    IF  X=2  AND  K=2  THEN 

1560      PRINT  "  SQUARE  ARRAY,  APEX  MOUNTED" 

1570  I  The  'L'loop  does  only  the  ape.-,  mounted  square 

1580     FOR  L=l  TO  N2 

1590        Rl=Real(L,L) 

1600        6=Itnag(L,L) 

1G10        Q  =  SQR(R1*R1+G*G  ) 

1620       IF  L'>PmaA  THEN  GOTO  1650 

1630       Amp(L)=Q/MaA 

1640       Ang(L>  =  ASN( ( (L-1  )  #68 . 6#R  )/ ( D*N  )  )  • 1 S0/P I 

1650  •    PRINT  L  ,Q/Ma,<  ,ASN(  (  (L-1  )*G8.6*R  )/(D*N  )  )#180/3.  1415927 

1660      NEXT  L 

1570  CALL  PolarpIt( Amp( ♦  ) ,Ang( ♦  )  ,PmaA  ) 

1680     IF  X-2  THEN  GOTO  2590 

1690  END  IF 

1700  I 

1710  I 

1720  FOR  1=1  TO  Bl 

1730     FOR  J=l  TO  B2 

1740     IF  Bl=l  THEN  P=J 

1750     IF  B2=l  THEN  P=I 

1760      Ri=Re3l(I,J) 

1770      6=Imag( I  ,J  ) 

1730      Q  =  SQR(R1*R1+G»6  ) 

1790       IF  P-'PmaA  THEN  GOTO  1840 

1800      Amp(P)=Q/Max 

1810      Ang(P  )  =  ASN( ( (P-1  )»68 . 6 'R  )/ ( D*N  )  )• 1 80/P I 

1820  I    PRINT  P  ,Q/Max  ,A5N( ( (P-1  )*68.6*R)/(D»N)  )*180/3. 1415927 

1830     NEXT  J 

1840  NEXT  I 

1850  I 

1860  CALL  Polarplt( Amp( *  )  ,Ang( •  )  ,Pmax  ) 

1870  I 

1880  B2=B1 

1890  Bl=l 

1 900  NEXT  K 

1910  ! 

1920  IF  A=3  THEN  GOTO  2590 

1930  I 

1940  I  This  section  centers  the  FFT  onto  the  output  grid. 

1950  1 

1960  IF  C=7  THEN   !  USE  ONLY  TO  GET  ONE  QUADRANT  OF  THE  GRID 

1970     GCLEAR 

1980     GINIT 
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1  99^3 

2010 
2020 
2030 
2040 
2050 
20G0 
2070 
2080 
2090 
2100 
2110 
2120 
2  I  30 
2140 
2150 
21B0 
2170 
2180 
21  90 
2200 
2210 
2220 


GRhPHICS  ON 

^JIEUPORT  30  ,  130  ,0  ,  100 

WINDOW  1  ,N2  ,  1  ,N2 

B!=N2-1 

B2  =  0 

GOTO  2340 
END  IF 
I 

I  This  section  centers  the  output 
I 

FOR  1=1  TO  N2-1 
FOR  J=l  TO  N2-1 

RI=Real( I  ,J  ) 

Id  =  Imag( I  ,  J  ) 

Real<N2  +  I-1  ,N2  +  J-1  )  =  RI 

RedI(N2-I  ,N2  +  J-1  )  =  R1 

ReaI(N2-I  ,N2-J  )  =  RI 

Real(N2  +  I-l  ,N2-J  )  =  R1 

Imag(N2-I  ,N2-J  )  =  Ia 

Imag(N2  +  I-l  ,N2-J  )  =  Ia 

Inag(N2  +  I-l  ,N2  +  J-1  )  =  Ia 

Inag(N2-I  ,N2  +  J-1  )=Ia 
NEXT  J 
NEXT  I 


2240 
2250 
2260 
2270 
2280 
2290 
2500 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 
2390 
2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 


Bl-B 
B2  =  B 
I 

GCLEAR 
G  I  rj  I  T 

GRAPHICS  ON 
I 

VIEWPORT  30  ,  130  ,0  ,100 
WINDOW  N2-B  ,N2  +  B  .N2-B  ,N2  +  B 
I 

FOR  'l=N2-BI  TO  N2+B2 
FOR  J=N2-B1  TO  N2+B2 
Rl=ReaI(  I  ,J  ) 
Ia  =  Imag(  I  ,  J  ) 
Q  =  SQR(R1*R1  +  Ia»la  ) 
IF  A=1  THEN 

Z  =  (  12»Q+1  )/lG   i  USE  ONLY  FOR  PLOTTING  INPUT  POINTS 
END  IF 
IF  A=2  THEN 

Z  =  l+LGT(Q/MaA  )/D   i  0NLY  FOR  DECADE  INTENSITY  SPREAD 
END  IF 

IF  Zv0  THEN  Z=0 
IF  Z'l  THEN  Z=l 
AREA  INTENSITY  Z  ,Z  ,Z 
m^E    I  ,J 
RECTANGLE  1  ,1  .FILL 


65 


2500 

NEXT  J 

2510 

NEXT  I 

2520 

1 

2530 

PRINT  "The  5ize  of  the  output  uindouj  15",B,"X",B 

2540 

PRINT  "The  size  of  the  input  u/indow  15",N,"X",N 

2550 

PRINT  "The  radius  of  each  element  i5  " ,R 

2560 

1 

2570 

1  Pause  statements  allow  insertion  of  parameters  u/o 

2580 

'  redoing  the  FFT  calculations. 

2590 

PRINTER  IS  CRT 

2B00 

PAUSE 

2610 

PAUSE 

2620 

END 

2630 

1 

2640 

1  This  section  is  the  FFT  subroutine  (from  INFOTEK). 

2650 

1 

2660 

SUB  Fft( INTEGER  M) 

2670 

con  Re(  1:512)  ,Im(  1:512) 

2680 

1  FFT  SUBROUTINE  FROM  INFOTEK  29  AUG.  1985: 

2690 

1  USE  COMMON  FOR  FASTER  EXECUTION   5.09  sec  vs  3.59 

2700 

INTEGER  I  ,J  ,K  ,L  ,N,Ip  .Le.Lel 

2710 

N  =  2"M 

2720 

J  =  t 

2730 

FOR  1=1  TO  N-1 

2740 

IF  I<.J  THEN 

2750 

T  =  Re( J  ) 

2760 

Re( J  )  =  Re(  I  ) 

2770 

Re( I  )  =  T 

2780 

T=Im( J  > 

2790 

Im( J  )  =  Im(  I  ) 

2800 

Im(I)=T 

2810 

END  IF 

2820 

K=N  DIU  2 

2830 

WHILE  K<J 

2840 

J  =  J-K 

2850 

K=K  DIU  2 

2860 

END  WHILE 

2870 

J  =  J  +  K 

2880 

NEXT  I 

2890 

Le=1 

2900 

FOR  L=1  TO  fi 

2910 

Le1=Le 

2920 

Le=Le+Le 

2930 

Ure=t 

2940 

Uim=0 

2950 

Ang=PI/Le1 

2960 

Wre  =  COS(Ang  ) 

2970 

Wim  =  SIN( Ang  ) 

2980 

FOR  J=l  TO  Lei 

2990 

FOR  I=J  TO  N  STEP  Le 

3000 

Ip=I+Le1 
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3010  Tre  =  Re(  Ip )*Ure-Im(  Ip )»Uim 

3020  Tim  =  Re(  Ip  )*Uin+Im(  Ip  )*Ure 

5030  'Re(  Ip  )  =  Re(  I  )-Tre 

3040  In(  Ip  )  =  Im(  I  )-Tim 

3050  Re(  I  )  =  Re(  I  )  +  Tre 

3060  In(  I  )=In(  I  )  +  Tin 

3070         NEXT  I 

3030         T=Ure*Ure-Uim»Uin 

3090         Uim=Ure*WiPi  +  Uim*Wre 

3100        Ure=T 

3110       NEXT  J 

3120     NEXT  L 

3130     Np5Cdle=l ./N 

31  43  MriT    Re=    (Npscaie  )*Re 

3150     MAT  In=  (Npscaie  )*Ipi 

3160  SUBEND 

3170  I 

3180  I  This  section  is  the  polar  plot  subroutine  written 

3190  '  by  Paul  Davison. 

3200  I 

3210  SUB  Polarpit  (  Amp(  •  )  ,Ang(  »  )  ,Pma,v  ) 

3220  REM    ++++++++++++++++++++++++++++++++++++++++++++ 

3230  REM    THIS  MODULE  MAPS  THE  CALCULATED,  NORMALIZED 

3240  REM    D5C5  ARRAY  DS(*)  ON  LOG  NORMAL  POLAR  PLOTS. 

3250  REM    ++++++++++++++++++++++++++++++++++++++++++++ 

32B0  I 

3270  I  ++++  INITIALIZING  THE  PLOTTER  ++++++++++++++++++ 

3280  ! 

3290  GINIT 

3300  PLOTTER  IS  CRT  ,"  INTERNAL " 

3310  GRAPHICS  ON 

3320  GCLEAR 

3330  I  +  +  +  +  SCALING  THE  GRAPH  +  +  +  +  +  -»-  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  + 

3340  ! 

3350  UIEUPORT  25,100.20,95 

3360  SHOW  -40  ,40  ,-40  ,40 

3370  ! 

3380  LINE  TYPE  4 

3390  Ring=10 

3400  AXES  Ring,Ring,0  ,0  ,5  .5.2 

3410  LINE  TYPE  3 

3420  DE6 

3430  LORG  7 

3440  CSIZE  3 

3450  FOR  R=0  TO  4 

3460     LINE  TYPE  1 

3470     MOUE  0  ,R*Ring 

3480     LABEL  R«Ring-40 

3490     LINE  TYPE  3 

3500     FOR  Angle=0  TO  360 

3510       PLOT  R*Ring*COS( Angle  )  ,R*Ring*SIN( Angle  ) 
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3520 

NEXT  Angle 

3530 

PENUP 

3540 

NEXT  R 

3550 

FOR  Angle=30  TO  150  STEP  30 

35S0 

IF  Angle<>90  THEN 

3570 

PLOT  40 'COS (Angle  )  .40*5 IN( Angle ) ,-2 

3580 

DRAW  40 ♦COS (  180  +  AngIe  )  ,40*SIN(  180fAngle  ) 

3590 

END  IF 

3600 

NExT  Angle 

3610 

PENUP 

3520 

1 

3530 

LINE  TYPE  1 

3640 

DEG 

3650 

FOR  P=Pna^    TO  1  STEP  -1 

3B60 

R=40+20*LGT( Amp(P ) ) 

3670 

IF  R>=0  THEN 

3680 

PLOT  R*COS(Ang(P  )  )  .R«SIN(Ang(P  )  ) 

3690 

END  IF 

5700 

NEXT  P 

3710 

FOR  P=l  TO  Pmax 

3720 

R=40+20*LGT< Anp(P ) ) 

3730 

IF  R>=0  THEN 

3740 

PLOT  R*COS(Ang(P)  )  ,-R»SIN( Ang( P  )  ) 

3750 

END  IF 

3760 

NEXT  P 

3770 

1 

3780 

CLIP  OFF 

3790 

1 

3800 

LORG  2 

3810 

MOUE  40,0 

3820 

LABEL  "X" 

3830 

LORG  8 

3840 

MQiJE  -40  ,0 

3650 

LABEL  "-X  " 

3860 

LORG  6 

3870 

MO^E  0,-40 

3880 

LABEL  "-Y" 

3890 

1 

3900 

CLIP  ON 

3910 

! 

3920 

DUMP  GRAPHICS  t701 

3930 

GRAPHICS  OFF 

3940 

PRINT  CHR$(  12  ) 

3950 

SUBEND 
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